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SYNOPSIS 


Decision Making is the process of selecting an 
alternative among the available ones. When there arc 
more than one criteria for the selection, the problem is 
called multi-criteria problem. Most real life problems 
fall under this category. There are various methods to 
solve this class of problems. All the- so methods can be 
classified into four major catagories depending on the 
preference information provided by the decision maker end 
the stage at which it is provided. 

In this thesis we develop an algorithm for bicriterion 

posterior 

linear programs. The algorithm requires/ articulation of 
information. The algorithm generates all efficient extreme 
points in the objective space. Almost all the existing 
algorithms for these class of problems generate efficient 
points in the decision spoce. However, the algorithm proposed 
by Aneja and Hair [1] concentrates on the generation of 
efficient points in the objective spoce. Hence the proposed 
algorithm is compered with algorithm proposed by Aneja and 
Hair. Por both the algorithms computer codes in PORTRAIT 
17 were developed and implemented on DEC 1090, Computer 
system. 
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For computation'--.! simplicity, in.lo.nt; ; d transportation 
problems of various sizes --cro develops!. These were solved 
as linear programs. The size of the linear program varied 
from problems with 7 constraints 16 variable to la constraints 
49 variables. 

The execution times for the generation of one extreme 
efficient point for both the algorithms were compared for 
24 randomly generated problems of 4 different sizes. It was 
observed that the proposed algorithms is comput ation.all3 r highly 
efficient as compared to Anoja and Nair's algorithm. For the 
smallest size problem, i.e., 4x4 balanced transportation 
problem, the average execution time (average biased on 
6 problems) per extreme efficient point w as l/6 that of the 
Aneja and Hair's algorithm. For biggei sized problems, 
the computational efficiency was even still better. 


/ 



CHAPTER I 


IITT RO DF CTI OH 

Excision making is a process of selecting a possible 
course of action among th~ various available ones. There 
may be one or more crit .ria which wall hip justify the 
selection. If thoT'. arc more then one criterion thc-n the 
process is celled Multipl- Criteria Decision Making (MCDM). 
Formally the criteria are non-commensurable and conflicting. 
Of late there lias boon considerable research activity in 
th- area of multi-criteria d> nisi on making sinco most of the 
real life problems fall under this category. In the MCDM 
problems, the best alternative is to be selected keeping 
in the co no id ora t ion, the various constraints and criteria,. 

In the process of choosing the best alternative, the d 'cision 
maker is raouirod to give some preference information. Thn 
stage at which it is given is the basis on which the MCDM 
problems can be classified. Th classification of this 
problems is given in Table 1.1. 

The methods for whies, no articulation of pr^f er>.-nce 
information is given, assunr that the decision maker will 
be able- to accept the given solution. This assumption is 
not well justified. The methods whore the decision maker 
is required to give his preference information aprion, 
suffer from the disadvantage -that the decision maker gives 
this in an information vacuum. The methods which are based 
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on ordinal prcft-rencc information given apriori depend 
heavily on this and tho moment the decision maker finds that 
the information given at an earlier time does not hold any 
more a tho problem is to bo solved afresh using the new 
information. 

In the third method, which is classified as interactive 
method, the decision mak-,*r is part of the solution process. 
Thor, is no nc-.d of a priori preference inform- -ion arti- 
culation. The major advantage of this is that decision 
maker learns about the problem as ho proceeds. Since decision 
maker is part of the solution process, the prospects of the 
obtained solution being implemented, is bright. But the 
solution depends upon the accuracy of the- local preference 
the decision maker can indicate. This is one of hie disadvan- 
tages, these methods suffer from. Another disadvantage is 
that much more effort is required on the part of tho decision 
maker than other methods. 

The last category of methods, listed in Table 1.1 are 
the methods for which a posteriori artifulation of prefe- 
rence information is given. These; methods h^ve an advantage 
over other methods, because no assumption regarding tho 
decision maker’s utility function is required. 

The present thesis deals with this class of multi- 
criteria problems. Solving thu multi-criteria problem 
requires the generation o f all non— dominated (efficient) 

points. 
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The first generalization of the multi-criteria 
problem is the bicrit orion problems. If all the cons- 
traints, and both the criteria arc linear, than the 
problem is called bicriteria linear program. The 
importance of the bicritcrion problems can knot be over 
emphasized. Many situations could be formulated as bi cri- 
teria linear programs. For instance in a production 
system ono would like to minimize the cost of production 
and maximize the production. 

A survey of the literature indicates that most of 
the investigators [2,4,9,10] have approached the bicriteria 
problem through developing algorithms for generating all 
th" efficient points in the d-cision space. These algori- 
thms involve the consideration of n dimensional space, 
where n is the number of decision variables. However, 
the dimensionality of the solution space can be considerably 
reduced if instead of the decision space, the objective 
space is considered. In the present thesis, an algorithm 
which generates the extreme efficient points in the objective 
space is developed. The performance of this algorithm 
is compared with Anejn and Hair's algorithm [l] the 
only other algorithm reported in the literature which 
exploits objective space rather than oocision space for 
the generation of extreme efficient points. A s-t of 24 
randomly generated balanced transportation problems of 
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different sizes are considered for comparing the two algo- 
rithms. For "bo th tho algorithms , these transportation 
problems are structured as linear programmes. 

The thesis is organized in four chapters. Chapter II 
deals with problem statement, its formulation and the 
relevant literature survey. The proposed algorithm is 
presented in Chapter III. In Chapter IV, a comparison 
of the proposed algorithm with Anoia and Fair's algorithm 
is presented. Listing of cG.l the compute?* programmes 
developed for this work are given in Appendices A, B and C. 



CHAPTER II 


P FOB LEM STATEMENT AND LITERATURE SURVEY 


2.1 Problem Statement s 

Bicriterion Linear Programs (BLP ) may be stated 
as follows; 

Consider a Vector Minimum Problem (VMP) 

jjcx, Dx 

s.t. Ax ^b 

x 7/0 t 1 ) 

where x , C and D are n-dimensional vectors, A is a (m x n) 
matrix and b is m-dimensional vector. C and D are coefficient 
of the two criterion functions. A represents the technolo- 
gical coefficients of the activities and b is the require- 
ment vector. Ax -c b represents the constraints of the 
linear programme (LP). This statement is made without loss 
of generality, as maximization problems and greater than 
or equal to constraints can be readily incorporated using 
standard techniques of LP. 

let us denote the set of feasible solutions by 
s and its mapping by the two objectives into objective or 
criterion space by Y. Y is called the pay off set and 

is defined as 1 

y = s 2 ) j s x = Cot, z 2 = lte for some ifs] 
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A solution is called superior if it minimizes 
both the objectives simultaneously. In most real life 
situations superior solutions do not exist. The objectives 
are non-commensurable. Under those assumptions the 
decision maker is interested in efficient solutions or 
non-dominated solutions. For clarity we define efficient 
solution before proceeding further. 

Definition ; 

A feasible solution x° is said to be efficient 
(non— dominated ) if Gx i-C Cx° for some other feasible x 
implies that Dx°< Dx and vice versa. 

Essentially it means that no other feasible solu- 
tion is better with respect to both the criterion. Also 
by implication x £ S is dominated if there exists atleast 
one feasible solution x° such that Cx° Cx and Dx°< Dx. 
Thus, it is always wise to restrict ourselves to non- 
dominated solution. The solution to (l) can be interpreted 
as determination of the s.t of entire efficient points. 

2.2 Survey of Literature;. 

Most of the research in this area lias been concen- 
trated on getting all the efficient points in the decision 
space. One of the early papers in this area is due to 
Geoff rion [4]. He proposed a scalar maximization approach 
for obtaining all the efficient solutions. He converted 
the original problem into the following proolem, 

t 



8 


Problem P a ; 

Min. aCx + (l-a) Dx 

S . t . Ax 

x ]>0 (2) 

where a <£ (0, 1). The problem P a gives all efficient 
solutions if a is varied in the given range. He proves 
the necessary and sufficient conditions in the form of 
a theorem. Further, G-eoffrion suggests that his 

approach can be extended to multi objective programs with 
more than two criteria. Bacopoulas and Singer [2] adopted 
a constraint criteria approach to get all the efficient 
solution in the decision space. One criteria is converted 
into a constraint at certain level. By parametrically 
varying this level he proves that all efficient solution 
can be obtained. 

Steuer [9] has given .an algorithm which is a modi- 
fied version of the simplex algorithm. In the column selec- 
tion rule stage, the efficiency of the next solution is 
tested by introducing each possible candidate columns 
one by one. One LP is solved to test the efficiency or 
non domination of each point. They call this procedure as 
adjacent efficient basis algorithm. 

Yu and Zeleny [10] have given another modification 
to the simplex procedure. They r ef er their modified simplex 
procedure as adjacent non-dominated -basis approach. A 
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modified version of the column selection rule is used to 
examine non-domination. Essentially they use the c . of 

d 

each objective to test non-dominance. 

Both Steuer and Zeleny and Yu concentrate on 
generation of efficient points in the decision space, 
which requires more computing and book keeping. Furthc-r, 
they have addressed to multi criteria problem in general. 

Sndagopan [S] has suggested an interactive solution 
procedure for the bicritcria mathematical programs. 

Recently, Aneja. and Fair [1] have developed an algorithm 
basically for bicriterion transportation problems. They 
claim that their algorithm is applicable to bicriterion 
linear programs also. 

Their algorithm generates efficient extreme points 
in the objective space. The algorithm requires solving of 
(2K-3) linear programs if there are K ( >2) efficient 
extreme points in the objective space. As the present work 
considers the same problem and provides algorithm as com- 
pared to Aneja and Fair’s approach, the details of their 
algorithm are presented in Chapter III. 

The proposed solution methodology involves the 
generation of efficient extreme points in the criteria space, 
instead of decision space. This was done for the following 
reasons. Firstly, the dimensions of the solution space get 
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reduced if the objective space is considered. For instance 
in the case of bicriterion problems the objective space is 
two dimensional while the decision space is n-dimensional 
where n is the number of decision variables. Secondly, 
all extreme points in the criteria space correspond to 
one extreme point in the decision space. Moreover, there 
may be more than one extreme point in the decision space 
mapping onto the same extreme point in the objective space. 
Since the efficient frontier in the objective space is 
piecewise linear curve, it is enough to restrict ourselves 
to extreme efficient points. Thus, it is advantageous to 
look for efficient extreme points in the objective space. 



CHAPTER III 


GENERATION OR EFFICIENT POINTS 

3.1 Introduction ; 

In the multicriteria problems as already mentioned, 
one is interested in the generation of efficient points. 

It was also stated that efficient points in the objective 
space will be more useful. Because of the polyhedral 
property of the payoff set, the set of all efficient solu- 
tion is determined the moment the extreme points of the 
polyhedron are known. Since the primary motivation to 
this research is the algorithm proposed by Aneja and Nair[l], 
this will be discussed first, followed by the algorithm 
proposed by us. 

3.2 Aneja and Nair f s Algorithm ; 

3.2.1 Introduction ; 

Aneja and Nair proposed the following algorithm 
for generating the entire set of efficient points in the 
criteria space. Though the authors have concentrated on the 
determination of efficient extreme points for the bi-criteria 
transportation problem, they claim that their algorithm is 
equally applicable to bi-criteria IP problems. The algori- 
thm initially tests for the existence of a superior solution. 
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If such a solution does not exist, 'the linear programming 
problem is solved using objective functions which are 
positively weighted average of the objective functions. 
This process is carried out till all the efficient extreme 
points are determined. 


„ _ Z (D 
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3.2.2 Algorithm ; 

Step 0 ; Find z^ = Min {z ± j x£S) and z^ = Min (z 2 j 
and x Q S). Record (z^^, z^^and set k = 1. Similarly, 
find z5> 2 ^ = Min (z^ x£S) and z£ 2 ^ = Min (z 1 | z 2 = z^ 
and x^,S). If (z^, ) = (z^ 2 \ z)p ) , stop. Else 

record (z^, z^ 2 ^ ) and set k = k+1. Define sets I =|(l,2)j 
and E = 0 , and go to Step 1. 

Step 1 ; Choose an element (r, s) &- L and set 

I 

I 

let x be an optimal solution to the linear program. 

Min. 


a < r ’ s) = 

!4 S) 

-4 r) 

a (r» s) - 
& — 

|4 S) 

-4 r) 


^(r, s) Cz + J,r, s) Dx 


( 1 ) 


S.t Ax -4 b 
x ^ 0 

If there are alternate optimal choose an optimal 
solution x for which Cx is minimum. 

Let z- = Cx and z 2 = Dx 
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If (z , Zp) is equal either to (z| r ^, zi r ^) or 

( Z 1 * Z 2 S ) sot ® = E U|(r, s)( and go to Step 2. 

Otherwise record (zP^, zi k ^)such that zi k ^ = z 1 and 
( k ) *“ 

Z^ ±w = z 2 » and set k = k+1. 


L = iu|(r, k), (k, s)j 

Step 2 ; Set L = L - ^(r, s) j . 

goto Step 1. 


and goto Step 
If 1 = 0 stop. 


2 . 

Otherwise 


Aneja and Nair have proved the following results. 

Theorem 3.1 : A point z ^ = (z^\ z^) is an efficient 

extreme point in the objective space, if and only if 
(k) 

z is recorded by the algorithm. 

Theorem 3.2 s The algorithm is finite and reauires exactly 
(2k~3) iterations when it has kfe2) efficient extreme 
points. 

3.2.3 Numeric al Exampl e ; ' ‘ 

For better understanding of the Aneja and Nair's 
algorithm, a numerical example considered by them is dis- 
cussed. Table 3*1 contains the data for a balanced trans- 
portation problem with two cost criteria. The two numbers 
in each cell represent the cost coefficients. The number 
in the north west corner indicates c. .'s and " tlle 

X J 

in south-east comer denotes d^'s. 



Table 3.1- Data for the problem 


Supply 

8 



The solution is indicated in Table 3-2. Tho 
values given by the side denotes the objective function 
values. 

Table 3.2; Optimal solution with c. • cost 

J~ t J 

coeff ici ents 
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Similarly z = 167 (= z £ ( x Qi s ) . The corres- 
ponding solution is indicated in Table 3.3. The objective 
function values are indicated by the side. 

Table 3.3; Optimal solution with d^j costs 
coefficients „ 



(2 


(2 


208 

167 


Iteration 

No. 1 




Step 1; 

(r, 

s) & L = 

: (1, 2) 



4 r ’ 

= 98 

( = \ 265 

- I67i ) 


Jr, 

s) = 61 

(= (208 

- 143/ ) 

Min. 

4 r ’ 

s) Ox + 

(r, s) 

ci 2 

Dx 


S. t x £ s (2) 

The Table 3.4 indicates the modified cost 
coefficients S ^ c^ + a^ 1 *’ d^ for all j =1,... mn 

where m and n are the size of the transportation problem. 
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Table 3.4 s Cost for transportation problem 
in Iteration Wo. 1. 

Supply 
8 
19 
17 

Demand 11 3 14 16 



Table 3.5” Solution for problem corresponding 
to Table 3.4. 


v. 5; (3 


V 6 ) 


1 ) (16 


zp) = 156 


Solution to the problem (1) is shown in Table 3.5. 


156 , z 2 = Dx = 200 


Since (z_ , z 0 ) is not either equal to (zpP zp^) or 


. iiV/ \ -cj- | 5 r\ 

.(2) , W (2K 


(zp ; , (zp ; ), we have, (zp\ zp^) = (z^, Zg) . 
Wow, k = 3+1 = 4 and L = pi, 2), (1,3), (3,2)j. 


1 = L 


- 1(1, 2 )\ = p 


1,5) , (3 


, 2 ){. 


Iteration Wo. 2 


(r,~s) = (1, 3 ) 


Summary of further calculations are shown in the 


form of Table 3.6. 
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Table 3.6; Bicriteria Solution in the Objective 
Spac o. 



A computer program was developed to solve the 
problems by Aneja and lair's method. The solution shown 
in Table 3.6 required u’.vun linear program to 'be 
solved. 

In the following section, the proposed algorithm * 
for generating the efficient extreme points in the 
objective space is described. 

3 . 3 An Adjacent Extreme Point Tree -dur e; 

3.3.1 Background ; 

In Chapter II, it was pointed out that the deter- 
mination of efficient extreme points in objective space 
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rather than the decision space offers two important advan- 
tages in case of multicrit cria problems. Algorithm 
described in section 3.2 generates all efficient extreme 
points in the objective space. Even though the convexity 
of the efficient frontier is recognized by Aneja and Hair [1], 
it is not fully exploited. It was also mentioned that 2k -3 
linear programs need to be solved if k (k ^2) efficient 
extreme points exists. 

In this section, an algorithm is developed which 
attempts to secure an adjacent efficient extreme point in 
each pivot operation. As explained later, the algorithm 
may not always be successful in generating the adjacent 
efficient extreme points. However, this does not instruct 
the utility of the algorithm since the next extreme 
efficient point will be secured in few more pivot steps. 

Before wo describe the actual algorithm, certain prelimi- 
nary results wnich lay the foundation for the development 
of the algorithm are stated. 

3.3.2 Preliminary Results r 

Let us define the original problem a.s 

PO : Minimize ^OX, DX *J 

s.t. AX 4 b 

X ?0 

where C, and D are n-dimensional vectors representing the 
cost coefficients, of the two objectives. X is a 
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n-dimensional vector of variables. A is m x n matrix 
representing technological coefficients and b is m x 1 
vector of requirements. 

Consider the following single objective, mathe- 
matical programs 


PI : Minimise CX 
s.t. AX< b 
X > 0 


P2 : Minimise Di 

s.t. AX ^ b 

V 

X y 0 


let the optimal values of PI and P2 be and z 0 

± 2 

respectively. Problems P 2 and P - are defined as follows 


Minimise EX 
s.t. AX,Z b 
CX«r z n 


X^ 0 


Minimise CX 
s.t. AX ^ b 

IX h 
X *? 0 


let x* and x 2 oe the optimal solutions for the 

problems P z * and P z * respectively. Sadagopan [s ] has 
1 2 

suggested a procedure for reducing the PO problem to 

P 2 (P 2 ) problem with paranctrizntion of z 1 (z 2 ). The 
1 2 

procedure utilizes the following theorems. 


Theorem 3.5: 


In the optimal solution to the program P~ 
z 2 <£: [z*, IS* ] the constraint DX <^z 2 will be a 
ono i.e. if x solve P- then Dx = z 0 . 


where 

2 

binding 



Proof ; 

Assume the contrary i.e. 

Dx < z 2 -4 Di^ (3) 

2 * being the absolute minimum of C over S 

S 1 = < Cx ( 4 ) 

Two cases will bo considered. 

Ca se 1 ; let Cx^ - Cx (5) 

( 3 ) and ( 5 ) tat: on together contradict the fact that 

solves P * since x is a better solution. 

± z 1 

Case 2 % Lot CxJ 4 Cx ( 6 ) 

Consider the lino segment [x, x^] which lies in the 
convex set S. Over this line segment the linear functions 
C and D are non-increasing. Because of strict inequality 
in ( 6 ), there exists a feasible solution x r on the line 
segment [x, x*J such -that 

Dx < Dx* C z 0 

2 (7) 

Cx >Cx' > Cx* 

Inequality (7) indicates that x' Is a feasible solution 

to P- with Cx' <Cx which contradicts the optimality of x. 
z 2 


Hence Dx 
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Theorem 3.4 s 

x £ S is efficient if and only if x solves P— 

z 2 

where ? 2 G. [ z*, Px^ ]. 

The sufficiency and necessity of this theorem is 
discussed below. 

Sufficiency ; Let x solve P — (henceforth whenever not 

z 2 

mentioned z 2 [z 2 > Dx|^ ]). If x is to be efficient, then 
it needs to be shown that there does not exist any x € S 
such that, 

Case (i) Dx ^ Dx and Cx 4 Cx 

( 8 ) 

Case (ii) Cx Cx and Dx h Dx 

Case (i) ; Let there be an x (x£ S) such that Dx 4 Dx = z 2 . 

Solution x is a feasible solution to P- but not optimal since 

2 

the theorem 3.3 states that every optimal solution x° 
satisfies Dx° =** 2 . Hence Cx y Cx. Therefore, x is not 
efficient. 

Case (ii ) ; Consider some x Cir S such that Cx Cx, where 

x solves P- . Since x is optimal cannot be feasible. 

z 2 

(otherwise x would have solved P- ) . This results in 
Dx > Dx. Consequently x is not efficient. 

Thus x is efficient. 
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Moccssity t Let .5 Ids the Set of all efficient solutions. 
Suppose x <c E. Obviously Cx < Ox* us otherwise x* will 
dominate x. let Cx = z^. Asciion that x does not solve 
P— but some other :t ! £,3 solves P- . By theorem 3.3? 

Z _1 _ z i 

Oic = z 1 = Cx' . Sine 3 x' is optimal to P— , lx’ <" lx. 

1 __ Z 1 s 
However, because x belongs to B this implies that 


Dx’ y Dx. Hence Dx' 

v 


Dx, 


ConsGouontly x solves 


z 


1 


Thus it is possible to generate the entire s.-t of 

efficient solutions bv parametrically solving P— ; . 

z 2 

theorem 3.3, the generated solution will have specific 
levels of attainment (Zp), of the second criterion. 
Hadley [5] suggests that for the parametric variation 


of the right hand side of the constraint sot, dual °i~>plex 
algorithm can be used. 


3.3.3 Algorithm (Adjacent Extreme Point Procedure) : 

Step 0; Solve PI and P2. Lot ^ and Zp be the optimal 

solutions to PI and P? respectively. Solve P * and P * . 

1 Z 2 

let xj and x| represent the corresponding optimal 
solutions. 

Stop 1 ; let C and D b^ the ’elativc cost vectors of the 
objective function of problem PI whose cost vectors are C and 
D. (In the first iteration, if alternate optimum to PI exists 
then choose among the alternate optimum, that solution which 
minimizes the second objective with cost vector D). Determine 
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-in ad,]- cent extreme point solution in which die entering 
v- viable in Xj_ such that 


= Min 
j 



for c - 0 and d . > 0 


; 


(9) 


(In case of ti :• break thorn arbitrarily). 

Set k = k+l> record z[ k ^ = CX , z^ = DX 


If DX =-• Dx"^, stop. Otherwise j, go to r ^tep 1. 

Theorem 3.5 s 

Algorithm described above generates all efficient 
extreme points. 

Proof “ As mentioned in section 3.3.2, the basic idea of 
this algorithm is to generate all efficient extreme 
points by solving B- over Z 2 [ Z 2 J • O^y onfc 

constraint (Px< z 2 ) is added to PI. Hence it can be treated 
as right hand side parametrization problem for the added 
constraint. Dual simplex algorithm can be used to do the 
parametrization. 

Lot B be the- basis nvtrix associated with thi problem 
PI. The matrix A with objective function in cononical form 
is represented as follows: 


~°B °II 


V 

i 

Z I 

B N 


X N 


b I 

L j 


where B (m x m) and N (m x n-m) are basic and not basic 
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rnatric s of A. C is also divided into Cg (1 x n) and 
C IvI ^ x n “ m ) representing the basic r>nd non— basic cost 
vectors res pec fcively. With the addition of the constraint 
Dx a 2 , the constraints can be written as follows. 

= b 


BXg + -f 0 • s 


D B- T B + Vl'T + 1 * S = z 2 


(ID 


Let A ' and b' be denoted as in (12) and (13) respec— 


, . _ _* 

cively . B represents the basis of P~ . Then 

z 2 


A 


N 0 




( 12 ) 


b 


L Z 2 J 


(13) 


B 


B 0 


(14) 


*)r 

Taking tho inverse of B [6], we have 


B 


*-l 


B" 1 0 


-DgB -1 1 


(15) 


Thus 


# -1 A r is 


n$ follows 


*’** 


r~ 




- 

B" 1 0 


B 

N 0 


i 

B^N 0 

-Dgir 1 1 


Dg 

1 _ 


0 

W' 1 " 1 - 


( 16 ) 
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where Bg (ixm) end B^-(i x n-n) are "basic and non-basic cost 
vectors of D respectively. 


To determine the right hand side range of the added 
constraint, the system B* “ 1 b’ ^ 0 is to bo solved. It 
essentially corresponds to the condition 


B 


whore ^ 


D B B " 

w 


0 

1 


T 

i 

i 


OO 


£ 


0 


(17) 


-^1 


represents the range over which P— continues 

Z 2 


to represent the optimal basis. Hence tho range will be 

[ z.2’ z 2^’ w kere Zg is the upper limit and is the lower 
(k) 

limit of ^ by solving the system (17). 


The variable to leave the basis will be the slack 
variable corresponding to the last row and the variable to 
'nter the basis will be decided by the minimum ratio rule 
of the dual simplex method. Since B^-Bg B~^ F precisely 
corresponds to the relative cost associated with the 
objective function with cost .coefficients of D vector, the 
rule to decide the incoming variable in Step 1 of the 
proposed algorithm corresponds to tho minimum ratio rule. 


Theorem 5.6 : 

Algorithm described in Section 5.3.3 terminates in 
finite number of iterations. 
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-J TPP — * 111 simplex,* no basis is ever repeated. Xn 

fn.ct the previous basis becomes infeasible to successive 
ones. The number of extreme points in the objective space 

is finite. Hence the algorithm terminates in finite number 
of iterations. 

It should b-_ noted that there may exist tie for the 

entering variable at some iterations, i.e., more than 

ono variable has the minimum rotio of — c ./d ■ where c". 

, 3 0 3 

and d^ are dual costs of non-bnsic variables. The tie can 
be bromen arbitrarily. In such cases, next adjacent 
extreme efficient point may not be obtained in one pivot 
operation. In case degeneracy is. encountered at some 
iterations, the lexicographic rule of simplex method can 
be adopted to prevent cycling. 

3.3.4 Numerical Example ; 

The balanced transportation problem which was given 
as an illustrated example in Section 3.2.3, is discussed 
to bring out clearly, the features of the proposed algorithm. 
Data for the problem is given in Table 3.1. 

Solution ; 

Step 0 ; Solving PI and P2, we have z£ = 143, = 167. 

The corresponding solutions appear in Tabic 3.2 and Table 3.3. 
Iteration Ho. 1 ; 

Step 1 ; Starting with the solution given in Table 3.2, 
the dual costs c. . and d . . of the non-basic colls are 

X J X J 
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determined . They are shown in Table 3.7. 

Table 3.7; Dual costs for the transportation 
problems at iteration no. 1. 


1 

v . 
3 

2 
v . 
3 


C 5 ) 

( i) 

5 

-10 

• 

3 

-5 

© 

v 

3 

T 

-5 

(n 

\t-y 

5 

io 1 

5 

1 6 


c?' 

1 

2 

4 

4 




a 2 


0 0 


0 


2 -8 


( 1 ) 

1 

( 1 ) 

2 


143 

265 


4 


13 


9 


In Table 3.7> the numbers in the north west corns; 
indicates the dual cost c. i and the numbers in the south 
oast corner represents the dual cost d. ... u. and v. 
around the table represents and v. with respect to 


2 - 2 


c. . costs, uh and v . are u. and v. with respect to d . . 

13 1 3 1 3 J 

cost s. 


Cell (2,3) enters the basis in the- next iteration j 

since the ratio of - c. ./d . is minimum for this non-basic ' 

1 J 1 .j j 

cell. 

low k = 0+i , = 143, = 265. j 

Since is not equal to 167 (Dxp go to Step 1. j 

I 

! 

| 

i 

# | 

l 

I 
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Iteration No . 2 : 


Step 1 ; 
revised . 
function 


Now the dual costs of non— basic variables arc 
They appear in Table 3.8 with the objective 
values by the side. 


Table 3.8: Dual costs of non -bn sic cellgat Iteration No. 2. 



Further calculations are shown in Table 3.9. 


Table 3.9; Solution to the transportation problem using 
the proposed algorithm. 


Iteration 

No. 

/ 

Leaving 

Variable 

Entering 

Variable 

r 

SI. No. of 
extreme 
efficient 
point 

00 


1 

(2,4) 

(2,3) 

1 

(143, 269) 

2 

(1,1) 

(1,3) 

2 

(136, 200/ 

3 

(3,3) 

(3,2) 

3 

(176, 179) 

4 

(1,2) 

(2,2) 

4 i 

(186, 171) 

5 



5 i 

(208, 167) 
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The set of efficient extreme points were generated 
in 5 iterations. Only two line-r programs were solved 
and the rest of the calculations were only pivot operations. 
Whereas Aneja and lair method reauiros solving of seven 
linear programs to arrive °>t the same set of extreme 
efficient points. 


* 



CHAPTER 17 


RESULTS AND DISCUSSIONS 

Computer codes were developed for the proposed as 
well as the Aneja and Hair’s algorithms. The programmes were 
written in FORTRAN 17 .and implemented in DEC 1090 computer 
system. Since "both the algorithms require solving of 
linear programs, a standard linear program code from 
Internationa]. Mathematical and Statistical library was utilised. 
The listings of programmes for Aneja and Hair's algorithm 

f ’ 

and the proposed algorithm are given in Appendix A and S 
respectively. The execution tames for the randomly gene- 
rated problems of various sizes were recorded for both the . 

algorithms along with the number of extreme efficient points : 
generated. It needs to bo pointed out that the lexicographic ; 
rule required to prevent cycling due to degeneracy was not 
incorporated in the computer code. This results in repeti- 
tive generation of certain points. However, very few points j 

exhibited this characteristics. The tie for entering variable ! 

I 

in tho proposed algorithm is resolved using the strategy of 
steeper reduction. This means that of all tho non-bnsic j 

variables, that x } is chosen to enter the basis for which 
the product c 1 9 is maximum. Hero S represents the level at : 

J 

which x i cnicpB "tho "fao^sis* 

j 
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1 * 1 Generation of Random Problems ? 

In order to generate feasible o nd bounded linear 
programs, balanced transportation problems were generated. 

The cost coefficients generated were in the interval 1 to S. 

The demand and supply wore in the interval 1 to 99 for smaller 
problems and 1 to 15 for bigger problems. A separate Computer 
programme was developed for the generation of random problems. 
The programme listing is given in Appendix C. 


The randomly generated balanced transportation problem 
of size m x n, is converted into a linear program with 
m + n - 1 constraints and mn variables. For example, a 
4x4 balanced transportation problem results in a linear 
program with 7 cons train to and 16 variables. 

4 . P Computation Time s 


Computational experience for 24 randomly generated 
balanced transportation problems was gathered for both 
the algorithms. Problems of four different sizes viz., 
i. x 4, 5x5, 6x6 and 7x7 were generated. Further, six 
problems for each problem size wore consider'd. The compu- 
tational times for solving these problems (in case of both 
the algorithms) are given in Table 4.1. Por each problem 
size, Table \.2 gives the average execution times for the 
generation of one efficient point for both the algorithms. 
The l«st column of this table gives a ratio measure to 

j t pi f H 

compare the computational efficiency of g^two. .algor if 


te* Mow 


* 

# 


3^ iftfc 
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Table- 4.1; Execution times by the proposed Algorithm 
for problems examined. 


Problem 

Size 

No. of 
efficient 
extreme 
points 

Time taken 
for gene- 
ration of 
extreme effi- 
cient points 
by Anejn and 
Nair's algo. 

( t ime in 
10--3sec . ) 

llo. of points 
generated by 
the proposed 
algorithm. 

Time t°ken 
for genera- 
tion of 
efficient 
points 

( t ime in 
10-3soc.) 

4x4 

1 

8 

3397 

8 

394 


2 

2 

839 

2 

344 


3 

2 

2779 

2 

381 


4 

5 

2317 

6 

384 


5 

4 

2060 

4 

379 


6 

3 

1078 

3 

283 

5x5 

1 

2 

1645 

2 

66 5 

2 

4 

4008 

4 

770 


3 

5 

5260 

5 

810 


/i 

6 

7351 

6 

925 


5 

7 

8465 

7 

869 


6 

9 

9220 

9 

808 

6 x 6 

1 

4 

8418 

4 

1642 

2 

4 

8181 

5 

1781 


3 

8 

17112 

8 

I 860 


4 

9 

20723 

9 

2129 


5 

6 

10 

22685 

13 

2286 


13 

30853 

14 

2388 

7x7 

1 

6 

18155 

6 

2321 

2 

7 

28897 

7 

2337 


3 

4 

5 

6 

9 

9 

11 

30857 

29105 

38964 

11 

9 

12 

2588 

2326 

2686 


12 

4 2863 

14 

2957 
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for the generation of one efficient extreme point. Of the 
24 problems considered , for three problems, there was 
repetitive generation of the sane point at some iterations. 
For one problem of size 6x6 there were three repetitions 
while for two problems of 7 x 7 there were two repeta -ons 
each. As mentioned earlier, the lexicographic rule to avoid 
cycling due to degeneracy was not incorporated in the 
computer code and this r esul ted in the repetitive generation 
of these points. Further whenever a tie existed for the 
entering variable the next adjacent extreme efficient 
point was not generated in one pivot operation. In case of 
tic a point was obtained on the boundary line Joining 
the present extreme point and the next extreme adjacent 
efficient point. This phenomenon was observed in one 
problem of size 4 i 4, two problems size 6 x 6 and one 
problem of problem size 7x7. Only one boundary point was 
generated in each of these problems. 

The average execution tine for one efficient point 
is calculated by dividing the total elocution “tino for all 
the six problems of particular size by the number of effi- 
cient extreme points as determined by Aneja and Hair's 
-algorithm. The second and third column of Table 4.2 show 
the average execution tine for one efficient point by Aneja 
and Hair’s algorithm and the proposed algorithm. 



Table -1 . 2: Average execution tines for generation of 

one efficient point by both the algorithms . 


5To. of 
problems 

Size 

Average _t?liro 
for generation 
of one efficient 
point by Aneja 
and Fair's algo- 
rithm. 

(A) 

Average time 
for generation 
of one effieier> + 
point by propo- 
sed algorithm 

(B) 

I 

Ratio ; 

(A) /(B) 

i 

6 

4x4 

445.35 

77.32 

5.75 

! 6 

5x5 

1089.36 

146 . 88 

7.41 

6 

6x6 

22-19.38 

251.91 

8.92 

6 

7x7 

2777.48 

281.38 

9.87 






From Table 4.2, it is clear that the proposed algorithm 
is computationally superior to Aneja and Fair's algorithm. For 
the 4 x 4 problem, the later algorithm required approximately 
6 times more execution time. As the problem size increases, 
the proposed algorithm shows even better computational 
superiority over the Aneja and Fair’s algorithm. Anejo and 
Fair’s algorithm requires solving 2k-3 linear programs, where 
k ( ^2) efficient extreme point exists, whereas the proposed 
algorithm requires only two linear programs and a little 
more than k pivot operations. This explains the reason for 
the computational superiority of the proposed algorithm over 
the Aneja and Fair's algorithm. 
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It should be- noted that the above stated computer 

tines are based on solving transportation problems using LP 

code. Further computational time reduction can be 

envisaged if a good transportation code is utilized, 

instead of the genor.nl linen- programming code. This is 

due to the fact, the transportation algorithm will not 

require the calculation of tho dual cost for all the 

variables after each pivot operation. Only the u^ s 

and • 1 s along .dth tho corresponding row and column 

3 , 

invoxvod in tho pivot operation need to be computed 

again. j 

i 

4 . 3 Conclusions ; j 

I 

In this thesis, an algorithm which generates efficient ; 
points in the objective space is developed to solve bicri- j 

. i 

terion linear programs. [ 

f 

The algorithm exploits the convexity of the efficient 
frontier in the objective space. The original bi criterion 
problem is reduced to a problem of parametrization using 
the results of Sadagopan [8]. This resulted in the 
proposed algorithm to require solving of two linear programs ; 
and a little over k pivot operations where k is tho number 
of efficient points. On the other hand, Aneja and Fair's I 

algorithm needs the solution of 2k-3 linear programs as the , 

authors do not fully exploit the convexity of the efficient : 


frontier. 



The comparison of the average execution tine to 
generate one efficient extreme point by both the algorithms 
clearly indicated the computational superiority of the 
proposed algorithm over Ane.ja and hair's algorithm. For 
instance, the smallest problem i.e., 4x4 balanced trans- 
portation problem required only one sixth of the average 
execution time of Anuja and hair's algorithm. The computa- 
tional superiority is even more significant when larger 
problems are considered. 

For very large problems number of efficient points 
may be large. In these cases an interactive procedure could 
he developed, which incorporates the proposed algorithm, so 
that the search region of the efficient frontier is reduced. 

This algorithm as such can not be extended to problems 
with more than two criteria. However, there is a need to 
develop algorithms for generation of efficient points in 
the objective space for multi criteria problems. 
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APPENDIX A 


Program Listing for the Proposed Algorithm 

m m m m m m *0 p+ #w m m ** m mmm §#&***#** 


00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

0 B 00 

01300 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

02100 

02200 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

04100 

04200 

04300 

04400 

04500 

04600 

04700 

04800 

04900 

05000 

05100 

05200 

05300 

05400 

05500 

05600 

05700 

05800 

05900 

06000 

06100 

8S338 

06400 

06500 


C*** MAIN ROUTINE „ 

COMMON /Bl/M,N.MA(S0),MB(5O),ICOSTi(5O,5O).ICOST2(5O,5 
COMMON /B4/A(50,50),B(50),IA,RW(4000) ,IW(2000) 

COMMON /81 1 /BASIS ( 50 ) » I BASIS (50 ) 

COMMON /B13/IZCNT,IZ(50,2) 

COMMON /B6/XC50. 503 
COMMON /TO/INPUT, IOUT 
COMMON /IOl/IOUTl ,INPUT1 
COMMON /625/ALPHA ( 5 ) 

LOGICAL END 

INTEGER BASIS, DBASI5.TTME1 

DATA INPUT , IOUT , IQUTi , INPUT 1 /20 ,5,22*23/ 

DATA IA/50/ 

END*, FALSE. 

NINE=9 

OPEN (UN IT* INPUT* FILES' TP ANSP.DAT' ) 

OPEN (UNIfsiQUT, FILE® 'TRANSP, OUT') 


OPEN (UNITsiOUTi, FILE* 'TIMING. REP ' f ACCESS* APPEND ) 


C 

c 

c 

c 


401 


c ^ 

C*** 

c 


to 


ur iutm-xuyu# r i 3 ul- 

OPEN (UNIT=INPUT1 .FILE® TIMING. OLD *3 


1000 


TYPE 401 
FORMAT! 'Give the number 
ACCEPT * » NQPRQB 
DO 1200 IIJJ=1*N0PR0B 
CALL RTIME(TIMEl) 

CALL ZERO 

CALL DREAD(END) A 

IFCEND) GO TO 1000 

CALL SETUP 

CALL FIRST 

CALL SECOND 

TYPE 4.IZCNT 

CALL TiMECTIMEl ,IZCNT) 

CALL PRINTt NINE) 

NINES9 
GO TO 1 
STOP 

CALL UERTST 
CALL ZX1LP 
END 


of problems') 


FIRST *** 


C0MM0nBb?/M*M?MA(50) ,MB(50) ,ICQST1C 50,50 ),ICOST2( 50* « 
COMMON /B4/A?56,50)*BC5Q),IA,RWC40QO),IW(2000) 


COMMON /B6/XC50.50) 

COMMON /B13/IZCNT,IZ(50,2) 


DIMENSION PSOLC50) ,DS0L(50) ,C0ST(50) 

INTEGER ROW, COL 
MlsM+N-l?M2sM*N 
DO 10 I = 1 ,M2 
PSOL ( I ) =0 . 0 
DSOLCI)*0,0 
C0ST(I)*0.0 

ROW* (I*l)/N+t 
COLsI-C (ROW-l)*N) 

COSTC I )®"IC0ST1 (ROW, COL) 

S ALL ZX3 LP ( A * I A , B , COST *M2*0#M1,S, PSOL » DSOL * RW , I W , IER ) 
DO 15 I * 1*M 


DO 15 



A* 1 


06600 
06700 
06800 
06900 
07000 
07100 
07200 
07300 
07400 
07500 
07600 
07700 
07800 
07900 
08000 
08100 
08200 
08300 
08400 
08500 
08600 
08700 
08800 
08900 
09000 
09100 
09200 
09300 
09400 
09500 
09600 
09700 
>91 


15 


8I!8§ 


20 


CONTINUE 

B(M1+1)=«S 

CALL ZX3LP(A#IA,B,COST#M2,0,M1+1 ,S#PSOL,DSOL,RW, IW,IER 

CAM, MATFQM(PSOL) , % 

C ALL XCHANG C PSOL # COST ,X,M2, SUM ) 

CALL PRINTC2) 

TEMP1=0.0 
TEMP2s0,0 
DO 20 I = 1 # M 
DO 20 J = 1,N 

TEMP1=TEMP1+IC0ST1 (I.J)*X(I#J) 

TEMP2*TEMP2+ICOST2 ( I . J) *X (T # J) 

CONTINUE 

' " ' 5 


c 

c*** 

c 


IZ(1 ,n=TEMPl + 0, 
IZ(l,2)=TEMP2+0, 
IZCNT=1 

CALL PRINT (13) 

RETURN 

END 




SECOND 

SUBROUTINE SECOND 


10000 

10100 

10200 

10300 

10400 

10500 

10600 

10700 

10800 

10900 

11000 

11100 

11200 

11300 

11400 

11500 

11600 

11700 

11800 

1900 

2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3400 

3500 


100 


4 


505 

5 


10 


20 

30 



50,5 


/B13/IZCNT#IZC50,2) 
/B6/X(50#50) 


COMMON 

COMMON # dp/ a ww » ~> v # 

INTEGER ENTER #LEAVE 
LOGICAL STOP # UBCOND 
STOP*. FALSE. 

M1*M+N-UM2=M*N . % 

SmTBB58 T »3*l'f”” , mrnMn , 

CALL FDLEAV CENTER# LEAVE. UBCOND) 
IFCUBCOND) CALL ERRORCl), 

CALL PIVOTCMl #M2 # A, B» BASIS .ENTER #LEAVE) 
CALL COMPUT 
GO TO 5 
return 

END 


SUBROUTINE PIVOT CM #N, A #B. BASIS, KEYCOL.REYROW) 
DIMENSION A (50 #50) ,BC50) .BASIS (50) 


INTEGER BASIS „ . _ 

DENOMR® A( KEYROW , KEYCOL) 
DO 3 J = 1#N 
IFC ACKEYROW, 

DO 4 I =!.M 


,J).EQ.1.0)2,3 


IF(I ? EQ;k^YR0W)GO0„5 


ir vhvI.w) .NE.O.O) GO TO 3 
CONTINUE 
ICHANG=J 
GO TO 5 
CONTINUE 

FORMAT?’ SOME THING WRONG IN 51290 # ) 

A (KEYROW# 5 )*A (KEYROW # J ) /DENOMR 

BCKEYROw5=B(KEYROW)/DENOMR 

CONTINUE 

DO 30 I = 1 #M 

IF ( I, EO. KEYROW) GO TO 30 

k I ‘yilfejf?Kl?Cag).A ( KE«OW, J , 

S ( I ) is ( I ) - AC I , KEYCOL) *B (KEYROW) 
CONTINUE 
DO 35 I *1 #M 
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13600 

13700 

13800 

13900 

14000 

14100 

14200 

14300 

14400 

14500 

14600 

14700 

14800 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

15600 

15700 

15800 

15900 

16000 

16100 

16200 

16300 

16400 

16500 

16600 

16700 

16800 

16900 

17000 

17100 

i?188 

17400 

17500 

17600 

17700 

17800 

17900 

18000 

18100 

18200 

18300 

18400 

18500 

18600 

18700 

18800 

18900 

19000 

19100 

19200 

19300 

11188 

UK8 

19800 

19900 

20000 

20100 

20200 

20300 

20400 

20500 


l£(f .I qjKe£rQW5A(KEYRQW,KEYC0L)=1*0 

IF(BASI5(I) , EQ.ICHANG)BASIS(I )*KEYCOL 

CONTINUE 

RETURN 

END 


SUBROUTINE XCHANGCPSOL,COST,X,N,SUMI „„„„ 

COMMON /B1/NQSP,NUSHPT { MA(50) ,MBC50) ,IC0ST1 (50, 50) ,ICOS! 


cuhwjp* ¥ da / jvuor $ iiuDnr - 1 # na j # ‘‘jowu * $ a v v j 9 * 

DIMENSION PSOL(N) ,C0ST(N) ,XC50,50) 

INTEGER ROW, COL 

SUMsO.O 

O') 30 I = 1 , N 

R0W=(I-1)/N0SHRT +1 

COL=I-((RQW-1)*NOSHRT) 

X(RQW,CQL)=PSOL(I) 

SUMsSUM + PSOL (13 ♦COST (I 3 

CONTINUE 

RETURN 

END 

COMMON T /B1/M^N?HA(50),MBC50),ICOST1(50,50),ICOST2(50,50 
COMMON /B4/A(50,50),B(50),IA,RW(4000),IW(2000) 

COMMON /Bl 1 /BASIS (50) , IBASIS (50 j 
COMMON /B13/IZCNT,IZ(50,2) 

COMMON /B6/XC50 ,50) 

DO 10 I = 1,50 
MA(I)=0?MB(I)=0 
B(I)»0 
BASI5(I)=0 


IBASIS(I)=0 
IZCI,1)=0JIZ(I,2)=0 
DO 10 J = 1,50 


A(I,J)=0„Q 

x(i#a5*o,o 

to CONTINUE 

MsO|N«0?IZCNT=0 
DO 15 I * 1,4000 
RW( I ) *0,0 

15 CONTINUE 4 _ A/v 

DO 20 J a 1,2000 
20 IW(J)=0 

RETURN 
it m n 

SUBROUTINE TIME(TIME1 , IZCNT) ! 

COMMON /B25/ALPHA(5) 

COMMON /B1/m!njMA(50KMB(50),ICOST1(50,50),ICOST2(50,50| 
DIMENSION ALPHAO( 6) 

INTEGER TIME1 ,TIME2,TIMDTF : 

CALL RTIMECTIME2) | 

C 140 5bIt”4u|i;i|o (»LPHA0CI).Iji.6),IZCNT.TIMDIF 

150 FORMAT (6 AS , * I , ,1X,I7,1X,I7,2X) 

RETURN 

END 

C*** MATFOM , 




COMMON /B4/AC50,50),B(6w1,1A,kh|«uwu^,4.i 
COMMON /BU/BASIS(50) »g|&STS(50) 
DIMENSION BASE(50,50) ,PSOL(50) ,DS0I.(50) 


INTEGER BASIS, IPSOL (50), DBASIS 
M1«M+N-1 f M2=M*N 
DO 10 I » l » M2 

IPS0L(I)=PS0L(I)T.005 


JsO 

DO 20 I « 1 , M2 
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20600 

20700 

20800 

20900 

21000 

21100 

21200 

21300 

21400 

21500 

21600 

21700 

21800 

21900 

22000 

22100 

22200 

22300 

22400 

22500 

22700 

22800 

22900 


20 


25 


35 

40 

45 


C 

C 

c 

c 

c 


3000 C*** 

23100 C 
23200 
23300 
23400 
23500 
23600 

23700 4 

2 | 8°0 

23900 5 

24000 
24100 
24200 

24300 6 

24400 
24500 
24600 

24700 7 

24800 
24900 
25000 
25100 
25200 
25300 
25400 
25500 
25600 
25700 
25800 
25900 
26000 
26100 
26200 
26300 
26400 
26500 

urn 
26800 
26900 
27000 
27100 

!?§§§ 

275§E 


8 

o 


11 


12 

13 


IF C IPSOliCI ) «LE,0 « Q)GO TO 20 
J s J + 1 

BASISCJ)=I 

CONTINUE 

IFCJ.E0.MDGO TO 40 

DO 35 I = 1 , M2 

DO 25 J s DM1 

IFCBASISCJ).E0,I)GO TO 35 

CONTINUE 

BASISCMl 3=1 

GO TO 40 

CONTINUE 

DO 45 J = 1.M1 

D8ASISCJ)=BASIS(J) 

BASEC I * Jl^AC I , BASIS CJ) ) 

CALL INVERT (BASE, Ml) 

CALL M ATMLT ( BASE » A # B # M 1 ,M2) 

RETURN 

END 


SUBROUTINE INVERT (A#N) 

DIMENSION A(50,50).PIVR(50),PIVC(50),PIVEC50) .YC50) 

DO 4 I = 1.N 

PIVRC I) *0 

PIVC(I)=0 

K*1 

PIVECK)=0.0 
DO 9 I«DN 

IF (I^EQ.PIVR(L) ) GO TO 9 
CONTINUE 
DO 8 J=DN 
DO 7 L=1 

IF ( J.EQ.PIVC(L) ) GO TO 8 

CONTINUE 

AB = ABS ( A( I # J) ) 

CsABS(PIVECK)5 mrt n 
IF (AB.LT.C) GO TO 8 
T. I r I 
J J* J 

PIVE(K)sA(I,J) 

CONTINUE 

CONTINUE 

DO 13 I = D N 
B*ACI»JJ) __ __ 4 „ 

IF CI.EQ.ID GO TO 13 

64*1 , 

CONTINUE 

CONTINUE 

IF K (K,LE,N) GO TO 5 
DO 16 J=Lg 
DO 15 1=1 #N 
IRSPIVRCI) 

ic*pivc(I) 
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27600 
27700 
27800 
27900 
28000 
28100 
28200 
28300 
28400 
28500 
28600 
28700 
28800 
28900 
29000 
29100 
29200 
29300 
29400 
29500 
29600 
29700 
29800 
29900 
30000 
30100 
30200 
30300 
30400 
30500 
30600 
30700 
30800 
30900 
31000 
31100 
31200 
3 1300 
31400 
31500 
31600 
31700 
31800 
31900 
32000 
32100 
32200 
32300 
32400 
32500 
32600 
32700 
32800 
32900 
33000 

mn 

33300 
“3400 
13500 
3600 
13700 
33800 
33900 
34000 
34100 
34200 
34300 
34400 
34500 


15 

16 


19 

20 


C ** 
C*** 

c 


20 

25 

30 


35 

45 


55 


5 

15 


25 


Y(TC)=A(IR,J) 
DO 161=1, N 
A ( I , J)=Y ( 1) 

DO 20 1=1, N 

DO 19 J=1 , N 

JR=PIVR(J) 

JC=PIVC(J) 

Y(JR)=A(I,JC) 

DO 20 J=1,N 

ACT, J)=Y(J) 

RETURN 

END 


SUBROUTINE MATMLT (BASE, A ,B,M1,M2) fA , D4 

DIMENSION BASE (50 ,50) ,A(50,50) ,B(50) , A1 (50 ,50) ,B1 C50) 

DO 30 I = l, Ml 

DO 25 J = 1 ,M2 

SUM = 0,0 

DO 20 K = 1 , Ml 

SUM=SUM+BASE(I,K)*A(K,J) 

A1 (I , J)=SUM 

CONTINUE 

DO 45 I = 1 , Ml 


SUM=0 ,0 
DO 35 J = 

SUM = SUM 
81 ( I )=SUM 
DO 55 I = 
B(I)=BKI) 

DO 55 J = 1 ,M2 


1 , Ml 

+ BASE(I,J)*B(J) 


1 ,Ml 




R 
END 


COMMOlP/B4/A?5^?50) ,R(50)»IA»RWt4000)#IriC2000) 
COMMON /ail/BASIS(5O),IBA5IS(50) 

m A ** a w-i a a a- ft P % 1LB t A i IT n 


COMMON ’/Bl /M, N » MA (50) , MB (50), ICOSTK 50,50 ) #IC0ST2 (50,50 
COMMON /B6/XC50 , 50) 


COMMON /Bl3/IZCNT,IZf5Q,2) 
DIMENSION PS0L(50 5 , COST (50 ) 
INTEGER BAS IS, DBAS IS 
M1=M+N-UM2=M*N 


DO 5 I * l, M2 
CQST(I)*0.0 
PSOL( I )®0 ,0 

S8“iI n j = i.Mi 

C ALL ( XCH ANG (PSoE |[ COST , X , M2 , SUM ) 

irasmu o 

DO 25 I = 1 ,M 

suMilsuMi li c i , j ) ncpsT t ( I , J ) 

SUM2*SUM2+X(I, J)*ICOST2(I, J) 

IZCIZCNT, l)=SUMl+0,5 

CALL Z pRlfiT?2);CALL*PRINT(14) 

RETURN 

END 


:*** fndent *** 


SUBROUTINE FNDENT CENTER , STOP ) 

?OMMON 50? ,MB(50) ,IC0ST1 C50 ,50) , ICOST2 C50,5( 
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34600 

34700 

34800 

34900 

35000 

35100 

35200 

35300 

35400 

35500 

35600 

35700 

35800 

35900 

36000 

36100 

36222 

36300 

36400 

36500 

36600 

36700 

36800 

36900 

37000 

37100 

37200 

37300 

37400 

37500 

37600 

37700 

37800 

37900 

38000 

38100 

38200 

38300 

38400 

38500 

38600 

38700 

38800 

38900 

39000 

39100 

39200 

39300 

39400 

39500 

39600 

39700 

39800 

39900 

40000 

40100 

40200 

40300 

40400 

40500 

40600 

40700 

40800 

40900 

ttm 

41200 

41300 

41400 

41500 


15 


20 


C 

C 


22 

25 

28 


30 

35 

36 

40 

104 


40 

m 

103 


COMMON /B4/AC50,50),B( 50) ,IA,RW (4000), IW( 2000) 
COMMON /B11/BASISC50) .IBASISCSO) 

LOGICAL STOP 

DIMENSION DJBARC50) ,CJRAR (50) 

INTEGER ENTER, ROW, COL, BASIS, DBASTS 
STOPS, FALSE, 

M1=M+N-1?M2=M*N 

EMTERsO 

FMINS9999, 

DO 5 1=1 , M2 

CJBAR(I)=0,Q 

DJBAR(I)=0,0 

DO 40 0=1, M2 

DO 15 I = 1 , Ml 

IFCO|EOg6AS I Sa))GO TO 40 

DUMMY1=Q,0;DUMMY2=0.0 

DO 20 I =1 , Ml 

ROW= (BASIS ( I ) -1 ) /N+l 

C0L=BASIS(I)-C(R0W-1)*N) 

DUMMY1=DUMMY1+IC0ST1(R0W,C0L)*A(I,J) 

DUMMY2=DUMMY2+ICOST2 (ROW, COL) *A( I , J) 

CONTINUE 



ROW, COD- 


DJBAR(J)=ICQST2(R0W,C0L) 
IF (C JBAR( “ 


GO TO 22 


DUMMY 1 

DUMMY2 

(J).GE.O.AND.DJBAR(J) ,LT,0,0) 

IF (FMIN* ( -C J8ARC J) /DJBAR( 0) ))40,28,25 
ENTER=J , % 

58ISI3m;ifMW555«J,o 

DO 35 I = I, Ml 

IFCACI, ENTER), L£,0.0)GO TO 30 
RATI01*BCI)/A?I ,ENTER) 

IF(RATI01,GE,XMIN1)G0 TO 30 
XMINlsRATIOi 

IFCACI, J),LE.0,0)GO TO 35 
RATI01=B(I)/ACf ,0) - _ 

IF(RATiai.GE,XMIN2)G0 TO 35 
XMIN2=RATI01 

TFCXMIN1 ,LT, 9999,0. AND, XMIN2.LT. 9999.0) GO TO 36 
GO TO 40 

IF (CJBAR(J)*XMIN2.GE.DJBAR (ENTER) *XHIN1)MIN=J 

ENTERS MIN 
CONTINUE 

WRITE(I0UT,104)(BASIS(I) .1=1, Ml) 

FORMAT ( f BASIS! * #20(I5,1X) ) 

WRITE (I OUT, 101 ) (CJBARC J) ,0*1 »M2) 

WRITE (IOUT, 102) (DJBARCJ) , J*1 #M2) 

DO 40 I = 1 ,MJ 



C 

C*** 

c 


IF ( ENTER, NE,0) RETURN 
STOP*, TRUE, 

RETURN 

END 


*** 


CDLEAV 

SUBROUTINE FDLEAV ( ENTER , LEAVE , UBCOND) /«:« m 

COMMON /Bl/MjN x MAC5O) f MBC50),lCOSTlC5Oj5O)jICOST2(5O,5l 
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41600 

41700 

41800 

41900 

42000 

42100 

42200 

42300 

42400 

42500 

42600 

42700 

42800 

42900 

43000 

43100 

43200 

43300 

43400 

43500 

43600 

43700 

43800 

43900 

44000 

44100 

44200 

44300 

44400 

44500 

44600 

44700 

44800 

44900 

45000 

45100 

wt 

45400 

45500 

45600 

45700 

45800 

45900 

46000 

46100 

46200 

46300 

46400 

46500 

46600 

46700 

46800 

46900 

47000 

47100 

47200 

47300 

47400 

47500 

47600 

47700 

47800 

47900 

48000 

48100 

48200 

48300 

48400 

48500 


10 


c 

c*** 

c 


LOGICAL UBCQND 

XMI N=9999 , 0 } UBCOND= . TRUE. 

M1=M+N~1?M2=M*N 

DO 10 I = t, Ml 

IF CACD ENTER) .LE.0.03G0 TO 10 
RATI 0=13 ( D / A C I » ENTER ) 

IF (RATIO.GE.XMIN)GO TO 
XMIN=RATIO 
UBC0ND=. FALSE, 

LEAVE=I 
CONTINUE 
RETURN 
END 


10 


10 

101 


15 

102 

20 


C*** 

c 


SUBROUTINE ERROR (I) 

COMMON / 10/ INPUT, I OUT 
GO TO (10,15,10) I 

WRITECIOUT, 101) . 

F0RMATC15X, 'INFEASIBLE SOLUTION OR WRONG DATA*//) 
GO TO 20 

WRITE ( I OUT , 102 ) „ , 

FORMATC 15X, 'UNBOUNDED SOLUTION*//) 

X=-l 

ROOT=SORTCX) 

END 


SUBROUTINE DREAD *** 


112 


155 


156 

101 


SUBROUTINE DREAD(END) 

LOGICAL END 

COMMON /BI/MjnVmUcSO) ,MB( 50 ) , ICOSfl (50 ,50) , IC0ST2 (50 , 50 
COMMON /B25/ALPHA(5) 

READClftPUf?! *2, END=100) (ALPHA (I ), 1=1 ,5) , M.INTO,N 
FORMAT ( 5 A5 ,I4,A2,I4) 

READ( INPUT , 101 ) (MACI), 1=1, M) 

READ (INPUT, 101 ) (MBCJ) , J=1 ,N) 

DO 155 J = 1 ,M 

READ( INPUT, l6t ) (IC0ST1CJ,K) ,K=1,N) 

‘ DM. 


113 


DO 156 J = 

RE ADC INPUT, 101 ) ( IC0ST2C J ,K)»K=1,N) 
FORMATC40I3) 

WRITE ( I OUT ,113) (ALPHA (I ) , 1=1 ,5 ) ,M , INTO, N 
F0RMAT(//////5A5,I4, A2,I4///) 

DO 15 11=1,2 




m 


102 

103 


FORSMffti^OxlHfl^lfiPul^OSTllAfRjrSiS^/iOX^DESTN 


(ICOSTl (1 ,J),J»1,N) 
(IC05T2C1 ,J),J=1,N) 


[5) 

iX, 'SOURCE 
= 2 , M 


1 ' 9X,20I5) 


104 

10 

15 


2)&RITE(lbUT,lQ4)t,ClCQST2(I,JDJ=l,N) 


106 

107 


108 

100 


f6rmAtcs 

DO 10 I 
IF (I I. EC 

ttlllrtu 

CONTINUE 

WRITECIOUT, 106) (MACD, 1=1, M) 

FORMAT C 1 OX , ' SUPPLY ' 7X,20I5) 

FORMAT Cl OX, 'DEMAND' ,7X,20I5) 

WRITECIOUT, 107) (MB (J) ,J=1,N) 

FORMATC 1HD30X, 'Solutions by Anela and Nalr method'//) 

RETURN 

E«D=. TRUE, 

RETURN 

END 
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48600 

48700 

48800 

48900 

49000 

49100 

49200 

49300 

49400 

49500 

49600 

49700 

49800 

49900 

50000 

50100 

50200 

50300 

50400 

50500 

50600 

50700 

50800 

50900 


;!?c g 


C 

c*** 

c 


to 


15 


r 

£*** 

C 


. _ .00 

51200 

51300 

51400 

51500 

51600 

51700 

51800 

51900 

52000 

52100 

52400 

52500 

52600 

52700 

52800 

52900 

53000 

53100 

53200 

53300 

53400 

53500 

53600 

53700 

53800 

53900 

54000 

54100 

54200 

54300 

54400 

54500 

54600 

54700 

54800 

54900 

55000 

55100 

55200 

55300 

55400 

55500 


10 

110 


20 


25 

120 

125 


30 


35 

130 

135 


40 

140 


50 

150 


245 

250 

255 


256 


60 

160 


SUBROUTINE SETUP 

COMMON /B1/M,N.MA(50) ,MB<50 ) , IC0ST1 (50 ,50 ) , IC0ST2 (50,50 

COMMON /B4/AC50.150) .BC50) f IA,RW(4000) ,IW(2Q0O) ,LPCNT 

M1=M+N-1 


M2=M*N 
DO 10 I = l f M 
BCI) =MA ( I ) 

DO 10 J = 1 ,N 
IJ= (I~1)*N + 0 
A(I . 103 = 1 ,0 
DO 15 I s M+l ,M1 
BCD =MB( I-M) 

DO 15 J a 1 ,M 
IJ=CO-;l)*N + CI-M) 
ACI ,103 = 1 ,0 
RETURN 
END 


SUBROUTINE PRINTCCODE) 

Common /B5/izic§0)]iz2C50) ,izcnt,ilcso,2) ,ilcnt ' 

COMMON /B1 /M^N^MAC 50) ,MBC50 ) , I COST 1 (50,50) , ICQST2 (50 ,501 

E89R8S /B4/A(io?lSof * IX » RW ( 4000 )»IW( 2000) #LPC«T | 

COMMON /B14/MINCQC(2),MINC0D(2) 

INTEGER CODE 



’’TURN 

WRITE ( IOUT, 120) CX( 1 , J) ,0=1 ,N) 

DO 25 I = 2,M „ , ; 

WRITE ( IOUT, 125) (X(I, 03,0=1 ,N) 

CONTINUE - . , 4 

E0RMATC5X, 'The solution matrix follows*//5X,20CF6.i,2X) 
FOPMAT(5X,20(F6.1 ,2X3) 

WRITE (I OUT ,130) (X(1,0),0=1,N) 

WRITE (IOUT, {35) (X(I,0),0=1,N) 

FORMAT (lx , 'The shippina matrix f oliows*/5X,20(F7 , 1,2X). 
F0RMAT(5X,20(F7, 1,2X3) 

RETURN 


FORMAT^ iOX* 'This is not an efficient point*//) 
RETURN 


WRITE ( IOUT, 150) 
FORMATC50X, * THE 


'//) 


PROBLEM IS TERMINATED 

F0RMAT(20x'*No^of^LP| solved? * ,I6/20X, **~ 
WRITECI0UT,i50) are . %I5 , /40X#29( ^ 


FORMAT (30X, 'Points follow? */30x , — — 

1 ftf O ^ f jf 

WRITE(IOUT,256) (IZl (I) , IZ2CI) ,1=1 ,IZCNT) 

FORMAT? C29X,I5,5X, 15/3 5 

wIlTEC IOUT, 160 ) IZCMT ,1Z1(IZCNT) , IZCNT«I^2( IZCNT) 
FQRMAT(20X» 'The efficient points are ZC',13, , 1 ) 
15X, 'Z( * ,13, ,2)= ,15) 


*//30x* * 


* *.I5 
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55600 

55700 

55800 

55900 

56000 

56100 

56200 

56300 

56400 

56500 

56600 

56700 

56800 

56900 

57000 

57100 

57200 

57300 

57400 

57500 

57600 

57700 

57800 

57900 

58000 

58100 

58200 

58300 

58400 


70 

170 


800 

805 


810 

815 


820 

10005 


1000 


900 

1105 


C * 
c*** 
c 


RETURN 

£Sn?^9y T 'l 7 °} I ?l (1) ' I Z2Cl),lZl(2) # IZ2C2) 

fOnMATCISX.’initUi extreme points follow'/20X# 'ZC1 #1)* 

L»dF )s ,I5 ' ZC2.1) = \t5,' Z(2,2) = M5//) 

if# X In# Pi 

WR I TEC I OUT #8 05) 

FQRMATC50X # * Solutions by parametrically varying r.h.sid 
WRI TEC TOUT #120)CXC1»J)#J=1#N) 

DO 810 I X 2,M 

WRITE CTOUT#125)CXCI#J)»J=lfN) 

CONTINUE 

IF C IZCNT==2 ) GO TO 820 


WRITECTOUT,815) IZ1 CIZCNT) #IZ2(IZCNT) . , 

pnoMftm r pay - * HIN C Cl # J)*XC I , J)* ,F10 , 4 # 'CQRRE5PNG D(I,J)*' 


F0RMATC20X, 

1 ,5X#F10.4) 

RETURN 

WRITECIQUT, 10005) . 

FORMAT C SOX# ’The primal is Inf easblble *//50x , 'The proble 
2ted and next problem pursued # ) 

RETURN 

WRITECIQUT, 1105) 

RETURN 

WRITE Cl OUT# 1105) 

F0RMATC50X# ’ WHOLE RANGE COVERED' 

RETURN 
END 


) 
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APPENDIX B 

tl2%m feU&SS I2S ggUZS Algorithm 


00100 

00200 

00300 

00400 

00500 

00600 

00700 

00800 

00900 

01000 

01100 

01200 

01129 

01400 

01500 

01600 

01700 

01800 

01900 

02000 

olio8 

02300 

02400 

02500 

02600 

02700 

02800 

02900 

03000 

03100 

03200 

03300 

03400 

03500 

03600 

03700 

03800 

03900 

04000 

8IJ88 

04300 

04400 

8SI88 

04700 

04800 

04900 

05000 

05100 

05200 

os: 

05' 
05500 
05600 

81188 

05900 

06000 

81188 

06300 

06400 


C 

c 

c 

c 


c 

c 


PROGRAM MAIN ROUTINE 

**** VARABLE DECLARATION „ x 

COMMON /B1/M,N.MA(50),MBC50). ICQST1 (50 , 50 ) .ICOST2C50.50) 
COMMON /B4/A?50, 150), 6(50), IA,R*(4QQO),IW(2OO0),LPCNT 
COMMON /B5/IZ1C50) , IZ2 ( 50 ) , IZCNT , IL ( 5Q , 2 ) ,ILCNT 
COMMON /IQ/ INPUT, I OUT 
COMMON /B25/ALPHAC5) 

COMMON /I01/I0UT1 

INTEGER TIME1.TIME2,TIMDIF _ 

LOGICAL END, ILEMTY , INFCON, SUPRIR, STOP 
DATA INPUT, IOUT,IOUT1/20, 21 ,22/ 

I A»50 

0 PEN (UN I T= INPUT ,DEVICE=DSK ,FILE= ' TBANSP. DAT ' ) 

OPEN (UNITsIOUT , FILEs 'TRANSP .OUT' ) „„„„„ 

OPEN ( UN ITsIOUT, FILES ‘ TRANSP ..OUT' f ACCESS* ' APPEND * J 
OPEN (UNIT=IQUTl , FILE* f TIMING . REP * , ACCESS* APPEND ) 

1 CALL RTIME(TIMEl) 

CALL ZERO /c . Mn ^ 

CALL DREAD (END) 

IF(END) _ GO..TO__15O0 


C _ 

c*** 

c 


CALL FEACHK(INFCON) 

IFCINFCON) CALL ERROR (1 ) 

SET UP A MATRIX WITH EOUALITY CONSTRAINTS 

CALL SETUP 


C*** FIND INITIAL EXTREME POINTS 

L CALL I N IT AL (SUPRIR) 

IF (SUPRIR) 11,10 
11 CALL TIMECTIMEi.l) 

CALL PRINT(t) 

GO TO 1 


C ^ 

c*** 

C 100 


CHOOSE NEXT IL ELEMENT. IF IL IS EMPTY TAKE NEXT PROBLEM 

Ball ?hos!lHlemty.ipostn) 

IF(ILEMTY) go to 1000 


bit . X 

6 CALL STEP21 ( IPQSTN , SUM l , SUM2 , IR , I S , STOP ) 

IF(STOP)GO TO 1 

EJJ* IF thI Iodutior^is not^effeicient do no record 

c*** ELSE RECORD. DO THIS BY STEP22 

c m iEim^b sm2 - lR - ls) 

1000 

" ■ TF1VE®5 

1500 CLOSE (UN IT* INPUT) 

B- 1 



06500 
06600 
06700 
06800 
06900 
07000 
07100 
07200 
07300 
07400 
07500 
07600 
07700 
07800 
07900 
08000 
08100 
08200 
08300 
08400 
08500 
08600 
08700 
08800 
08900 
09000 
,00 
j 00 
09300 
09400 
09500 

mi 

09800 
09900 
0000 
0100 
10200 
0300 
0400 
0500 
0600 
0700 
0800 
0900 
1000 

tm 

1300 
11400 

mi 

11700 

11800 

B888 

12100 

12200 

BI88 

12500 

12600 

12700 

12800 

12900 

13000 


CLOSE (UNIT® IOUT) 
CALL UERTST 
CALL ZX1LP 
END 


SUBROUTINE ZERO 


COMMON /Bi/M,N,MA(50),MB(50) , ICQST1 (50 ,50) ,ICOST2(50,50) 
COMMON /B4/A(50,150),B(50),IA,RW(400p) , IW ( 2000) , LPCNT 
COMMON /B5/1Z1C50) ,IZ2(50),IZCNT,IL(50,2) ,ILCNT 
COMMON /B1 1 /BASIS (50 ) ,IBAS!S(50) 

COMMON /B13/IZCNTD.IZ(50,2) 

COMMON /B6/X(5Q,50) 

DO 10 I = 1,50 
MACI)=0,*MB(I) = 0 
B ( I ) =0 ^ 

BASIS (13=0 
IBASIS ( 13=0 
IZ1 (I)=0 

IZ(I*1)=0,*IZ(I,2) = 0 
DO 10 J = 1,50 
A(I , J)=Q,0 
X(I,J)«0,0 
10 CONTINUE 

M=0,*N=0,*IZCNT=0 
ILCNT=0 
IZCNTD=0 
LPCNT=0 . 

DO 15 I * 1,4000 
RN(I)=O#0 
15 CONTINUE 4 

DO 20 J » 1,2000 
20 IW(J)=0 

RETURN 
END 
C 

SUBROUTINE TIME (TIME!, IZCNT3 
COMMON /B25/ALPHA (53 

COMMON /Bl/Mj N^MA (50), MB (503, IC0ST1 (50,50), IC0ST2 (50 , 503 

SSt?E(IoU|Tni0)(k ? HMX),I = l,5),«,N,TIHDIF,IZCNT 
150 F0RMAT(5A5,2X, 12, 'X # ,I2,2X,I7,1X,I3,20X3 
RETURN 
END 

C*** SUBROUTINE DREAD *** 

C 


SUBROUTINE DREAD(END) 
LOGICAL END 



112 


fflos /I?/MjNjMiti 0 ).M 8 t 50 ),IC 0 SnC 50 , 50 ), 

COMMON /B25/ALPHA(5) 

lI»DtliSPUT?fl2 ( END»l(IO)(StPHfta),I=1.5), «,XNIO,N 
FORMATtSAS ,I4,X2 . I4> 

READ ( INPUT, 101 3 (MACI 3 , 1=1 ,M) 

RE ADC INPUT, 101 3 CMB( J) , J=1 ,N) 

DO 155 J = 1,M 


IC0ST2(50,50) 
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13500 

13600 

13700 

13800 

13900 

14000 

14100 

14200 

14300 

14400 

14500 

14600 

14700 

14800 

14900 

15000 

15100 

15200 

15300 

15400 

15500 

156O0 

15700 

15800 

15900 

16000 

16100 

16200 

16300 

6400 

S288 

6700 

6800 

6900 

7000 

7100 

7200 

7300 

7400 

7500 

7600 

7700 

7800 

7900 

8000 

8100 

8200 

8300 

8400 

8500 

,8600 

8700 

8800 

,8900 

.9000 

9100 

.9200 

.9300 

-9400 

19500 

19600 

19700 

19800 

19900 

20000 

20100 

20200 

20300 

20400 


155 

156 
101 

113 


102 

103 


104 

10 

15 

106 

107 


108 

100 


c 

c*** 

c 


10 


15 


c 

c*** 


c*** 

c 


READtINPUT,l01)CICnS*Pl(J,K) ,K=t,N) 

DO 156 J = l.M 

READ (INPUT, 101 )CTCOST2(J,iO ,K=l,N) 
FURMATC40I3) 

WRITE (I OUT, 113) ( ALPHA ( I ), 1=1*5), M , INTO , N 
F0RMATC//////5A5,I4,A2,I4///) 

DO 15 11=1,2 

WRITEC IOUT ,102) TI, (1,1=1 ,N) 



1 '9X,20I5) 


3 ^5X, 2015^ 

FQRMAT(5X, 'SOURCE 
DO 10 X ss 2 M 
IF(II.EQT2)^RITE(IOUT,10f)I,(lCpS|2CI,J) , J=1,N) 

IF( I I , EQ, 1 ) WRITE ( IOUT, 104)1 , C ICQ5T1 (I,J),J=1,N) 

FORMAT (10X,I4,9X,20I5/) 

CONTINUE 

CONTINUE 

WRITE ( IOUT, 106) (MA ( I) , 1=1 ,M) 

FQRMATC10X, 'SUPPLY ' 7X,2015) 

F0RMAT(10X, 'DEMAND ',7X, 2015) 

WRITEC IOUT, 107)(MB(J),J=1 ,N) 

WRITE(IOUT,108) , , , . . j * />. 

FORMAT ( 1H1 ,30X, 'Solutions by Aneja and Nair method'//) 
RETURN 
END=,TRUE, 

RETURN 

END 


DESTN „ «•> ' 


C0S«a» T /Bf/SS}C50),MB(50),ICOSTl(50,50).ICOS12(5O,50) 
COMMON /B4/A(50, 150) ,B(50) ,IA,RW(4Q00) ,1«(2000) , LPCNT 
MlsM+N-1 
M2=M*N 

DO 10 I = 1,M 
B(I)=MA( % ) 

DO 10 J = 1 ,N 
XJ= (I«l)*N + J 
A(I,IJ)3l.O 
DO 15 I = M+lfMl 
B(I)=MB(I*M) 

DO 15 3 * 1 ,M 
IJ*(J-l)fN i (I-M) 

A(I, IJ)=1 ,0 

RETURN 

END 


C0MM0N^/B1/m!!n3mA( 50)^MB^50 ) , ICOSTI (5Q ,50) ,IC0ST2 (50*50) 

Son /B4/A?s6^150?i^( 50 ) » I y??50 * 2?^ 

COMMON /B5/IZ1 (50) .IZ2C50) , IZCN1,ID(50,2) ,Il#CNT 

COMMON /B6/X(50.50) ncmri*;ni 

DIMENSION C0STC150) ,PSOL( 150) ,DSOlt(150) 

TWD RIR 


Ml*M+N-l 
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20500 

20600 

20700 

20800 

20900 

21000 

21100 

21200 

21300 

21400 

21500 

21600 

21700 

21800 


IS 


*000 
22100 
22200 
22300 
22400 
22500 
22600 
22700 
22800 
22900 
23000 
23100 
23200 
23300 
23400 
23500 
23600 
23700 
23800 
23900 
24000 
24100 
24200 
24300 
24400 
24500 
24600 
24700 
24800 
24900 
25000 
25100 
25200 
25300 
25400 
25500 
25600 
25700 
25800 
25900 
26000 
26100 
26200 
26300 
26400 
26500 
26600 
26700 
26800 
26900 
27000 

1 ?IB 8 

27300 

27400 


25 


28 


too 


c*** 

c 


M2=M*N 

DO 100 11=1,2 
DO 15 I =1,M 
DO 15 J =1,N 

ij=a-n*N+j 

IFCII.e 6.2)COSTCIJ)=-ICOST2(I,J) 
iF(ii.EQ.ncasT(ij)=-icosTi ci,j> 

CONTINUE 

CALL ZX3LPCA.IA,B,CUST,M2,0,M1 ,S,PSOL,DSOL,RW,TW,IER) 
LPCNT=LPCnT+1 

CALL XCHANG(PS0L»CnST,X,M2,T0TC0S) 

TW0=2 

CALL PRINTCTWO) 

IFC1I,EQ.1)IZ1(II)*-TOTCOS 

IF(II,EQ,2)IZ2(II)=-TOTCOS 

BCMl+i)=-TOTCOS 

DO 25 I = 1 , M 

DO 25 J s 1,N 

IJ=CI-1)*N +J 

A(M1+l.IJ) = IC0ST2a,J) 

IFCII.EQ. 1)ACM1+1,IJ)=IC0ST1(I,J) 

COST(IJ)=-ICOST2CI,J) 

CONTINUE* 25 C0STaj)= * IC0ST1(I ' J) 

LPCNT=LPCNT+1 

CALL ZX3LP(A,XA / B,C0ST,M2,0,Ml+l-S,PS0L,DS0L,RW f IW,IER) 
CALL XCHANG C PSQLf COST ,X,M2,T0TC0S) 

CONTINUE , . 

IF(II.EQ.1)IZ2(II)*-TQTC0S 
IFCII # EQ,2)IZ1CII)»*T0TC0S 
T«in 35 2 

CALL PRINTCTWO) 

CONTINUE 
ILCNT=1 
ILC1 , 1 )*1 
ILC1 , 2)*2 
IZCNT=2 

IFCIZl(iKEO^!zi(2).AND.XZ2Cl).EO,IZ2(2)) SUPRIR*.TRUE. 

RETURN 

END 


SUBROUTINE STEP21 ( IPOSTN , SUMt , SUM2 ,IR,IS, STOP) 

§§Ha§K ^ll^fI?M»PP4«I(^)jICoJll(io:loijlCOSI2C50.50) 

COMMON /B4/M 56 , ISO) ,6(50),I^,RW(4000),IWC^000) , LPCNT 
COMMON /B6/X(50,50) 

COMMON /B25/ALPHAC5). _ ncmMI! ., 

DIMENSION COST C 150) , PSOLC 150) #DSOLC150) 

LOGICAL STOP 
INTEGER THREE ,TWO 
STOP ® .FALSE * 

THREE * 3 
TW0®2 

IRsILCIPOSTN.l) 



10 


DO 10 I = 1 ,M 
DO 10 J * . 1#N 

COST(lil*IAl *IC0ST1 (I» J) + IA2*IC0ST2U,J) 
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27500 

27600 

27700 

27800 

27900 

28000 

28100 

28200 

28300 

28400 

28500 

28600 

28700 

28800 

28900 

29000 

29100 

29200 

29300 

29400 

29500 

29600 

21288 

29900 

30000 

30100 

30200 

30300 

30400 

30500 

30600 

30700 

30800 

30900 

31000 

31100 

31200 

31300 

31400 

31500 

31600 

31700 

31800 

31900 

32000 

11488 

32300 

32400 

m°o° 0 

32700 

32800 

32900 

33000 

33100 

33200 

33300 

33400 

33500 

33600 

11288 

33900 

34000 

11188 


15 

25 

35 

30 


50 

55 

601 


C*** 

c 


c 25 
C*** 
25 
20 


CQST(1)=-C0ST(I) 

CALL PRINT (THREE) 

M1=M+M"1 

M2=M*N 

LPCNT=LPCNT+1 

CALL ZX3LP(A,IA,B,COST,M2,0,M1 ,TOTCQS,PSOL,DSQL,RW,IW, JER) 
DO 35 I s 1 ,M 
DO 35 J = 1 * N 
IJs(I-l)*N+J 

A(M+N,IJ)=IA1*IC0ST1(I,J) + IA2*ICOST2CI, J) 

COST (IJ)=»IC0ST1(I,J) 

B(M+N)*-TOTCOS 

LPCNT=LPCNTti 

CALL ZX3LPCA,IA,B,COST,M2,0 # Mm,S,PSQL,DSOL,RW,IW,IER) 
TFUER,EQ.U3)G0 TO 55 
SUM1 =0 « Oj SUM2»0 ,0 

CALL XCHANG(PS0L,C0ST,X,M2,T0TC0S) 

CALL PRXNT(TWO) 

DO 50 X = t,M 
DO 50 J = 1 ,N 

SUM1=SUM1+IC0ST1(I,J)*X(I,J) 

SUM2*SUM2+IC05T2(I,J)*X(I,J) 

RETURN 
STOP=, TRUE. 

TYPE 60 1 # (ALPHA (III),IIIsl,5) ^ 

FQRMATC//5A5# ' IS INFEASIBLE*// ) 

CALL PRINTC14) 

RETURN 

END 


SUBROUTINE SETP22(SUM1,SUM2#XR»IS) 

COMMON /B5/IZ1 (50 ) . IZ2 ( 50 ) , IZCNT , XL (50 , 2) , ILCNT 
COMMON /B6/X(50,5Q) 

REAL SUM1 ,SUM2 
INTEGER FOUR # SIX 
FQUR=4?SIX=6 
ISUM1=SUM1 ? ISUM2=SUM2 v % 

DIFR1*ABS(SUM1-IZ1(XR)) 

DIFR2*ABS(SUM2-IZ2(|R)) 

OIFSlsABS(SUMl-IZlCIS)) 

DIFS2«ABS (SUM2"IZ2 ( IS) ) „ , . Mn nTK . co TT J 

IF ( (DXFR1 «LT. 1 .AND.DIFR2.LT.1 ) .OR . (DIP SI » LT» 1 .AND,DIFS2.I*T, 1 

1 ) )25 . 20 

C* hi r rRIMT ( FOUR 1 
WRITES THAT THIS POINT IS NOT EFFICIENT POINT 
RETURN 


C*** 


*** 


IZCNTsIZCNT+l w> 

IZ1CIZCNT)*ISUM1 

IZ2CIZCNT)sISUM2 

ILCNT*ILCNT+1 

ILCILCNT# 1 )*IR 

IL lILCNT » 2 5 sJZC NT 

ILCNT=ILCNT+1 

RE?&RDS C THE EFFICIENT POINT 
CALL PRINTC2) 

CALL PRINT(SIX) 

RETURN 

END 



34500 

34600 

34700 

34800 

34900 

10 

35000 

35100 

20 

35200 

35300 

35400 

35500 

C ^ 

35600 

c*** 

35700 

c 

35800 

35900 

36000 

36100 

36200 

36300 

36400 

10 

36500 

101 

36600 

367O0 

15 

36800 

102 

36900 

37000 

37100 

c 

37200 

c*** 

37300 

c 

37400 

37500 

37600 

37700 

37800 

37900 

3*000 

38100 

38200 

c 


.1300 
3*400 
3*500 
3*600 
38700 
38800 
38900 
39000 
39100 
39200 
39300 
39400 
39500 
39600 
39700 
39800 
39900 
40000 
40100 
40200 
40300 
40400 
40500 
40600 
40700 
40800 
40900 
41000 
41100 
41200 
41300 
41400 


10 

110 


20 


25 

120 

125 

30 


35 

130 

135 

40 

140 

50 

ISO 

245 

250 


SUBROUTINE SETP3 UR, IS) 

COMMON /B5/IZ1C50),IZ2(50),IZCNT,IL(5G,2) ,ILCNT 

IFC(ILUu5!eQUR).AND. CJLCI,2)«EQ,IS))GQ TO 20 

CONTINUE 

RETURN 

ILCI,1)=0 

IL ( I , 2 ) =0 

RETURN 

END 


SUBROUTINE ERROR(I) 

COMMON /10/INPUT, IOUT 
GO TO (10,15,10) I 

WRITE ( IOUT ,101) „ „ „ mm . # 

F0RMATC15X, UNFEASIBLE SOLUTION OR WRONG DATA'//) 
RETURN 

F0Rm!I(15X^UNB0UNDED SOLUTION'//) 

RETURN 

END 


SUBROUTINE PRINT(CODE) 

gg&Bgg /B5/IZl(§0)!!lZ2(50),IZCNT,Iti(5Q,2),I|,CNT 

COMMON /Bl/M?«“BA(50),MBt50),ICOSll(50,50),XCOST2(50,50) 

S8SS8S /B4/A( 50^lSof #B(|o ] , XX , RW (4000 ) *IW(2000) , LPCNT 

COMMON /B14/M1NC0C(2) ,M1NC0D(2) 

INTEGER CODE 

GO A TO(10?20;30?So?50 f 60, 70. 800, 900. 1000, 1000, 1000,1 000, 820)C 

WRITE!lo6T,ilo5;i£lU),fz2fn,fzic5) f IZ2(2) , 7(1 

FORM AT (15X, 'Superior solution exists /20X, ZCl.l)- 15, ZCl, 
1,15,' Z (2 , 1 )® , 15 , Z(2,2)=',I5//3 

WRITE U OUT, 120) CX ( 1 , J) , J*1 , » ) 

WRITE (f OUT, f 25) (XCI, J) , J-l ,N) 

CONTINUE . 

FQRMATC5X, 'The solu,--. 

F0RMAT(5X,20(F6,1 ,2X) ) 

WRITE (IOUT, 130) (XCI , J) , J»l ,N) 

WRITE (IOUT, 1 35) (XCI, J) , J=1 ,N) 

FORMAT (lx, 'The Shipping matrix follows ' /5X,20 (F7 , 1 , 2X) //) 
FORMAT (5xJ 20 (F7 , 1 , 2X) ) 

RETURN 

FORMATClOX^Thls is not an efficient point'//) 

RETURN _ 

FORMAT (§0X^ THE PROBLEM IS TERMINATED '//) 

solved* M6/20X, ~ //> 

fficient points are* ' , I5,/4QX,29C *„')// 
B- 6 


FORMAT (5X The solution matrix tol lows '/ /5X.20 (F6 • l » 2X) / ) 
FORMAT" *' 4 

RETURN 




41500 

41600 

41700 

41800 

41900 

42000 

42100 

42200 

42300 

42400 

42500 

42600 

42700 

42800 

42900 

43000 

43100 

43200 

III88 

43500 

43600 

43700 

43800 

43900 

44000 

44100 

44200 

44300 

44400 

44500 

44600 

44700 

44800 

44900 

45000 

45100 

5200 

45300 

45400 

45500 

45600 

45700 

45800 

45900 

46000 

46100 

46200 

46300 

46400 

46500 

46600 

46700 

46800 

46900 

47000 

47100 

47200 

47300 

47400 

47500 

47600 

47700 

47800 

47900 

48000 

48100 

48200 

48300 

48400 


255 


256 

60 

160 


7 < 
170 


800 

805 


810 

815 


820 

10005 


1000 

900 

1105 


*** 


10 

15 
C _ 

c*** 

c 


'//30x, 'One',8x, 'T 


,15, 


WRITE ( I OUT# 255 3 

FQRMATC30X, 'Points follows V30x, 

1 $$0 f / / ) 

WRITEC IOUT ,256)(IZ1(1),IZ2(I),I=1, IZCNT) 

FORMATC(29X,15,5X,T5/)) 

RETURN 

WRITE (IOUT, 160) IZCNT, IZ1 (IZCNT) , IZCNT, IZ2( IZCNT) 

FORM AT (2 OX, ’The efficient points are ZC%I3,%1 U ' 
lSX^'ZC ' ,I3» • ,2) = ',15) 

WRITE ( IOUT fl70)IZl(l),IZ2(l),IZiC2),IZ2C2) , ,, ^ 

FORMAT ( 15X i Initial extreme points follow* /20X, 'Z(l ,1)='I5, 

1 ' Z(1,2)=‘,I5,' Z(2»1)=*,I5,' ZC2,2)*',I5//) 

RETURN 

WRI TE ( IOUT , 805 ) 

FQRMATC50X, 'Solutions by parametrically varying r ,h„side'//) 
WRITE(IOUT,120) (XC1 # J),J=1,N) 

DO 810 I = 2 ,M 

WRITE( IOUT, 1 25 )(X(T,J),J=1,N) 

CONTINUE 

IFCIZCNT=»2)GO TO 820 v , „ 

WRITE ( IOUT ,815) IZKTZCNT) ,IZ2CIZCN|) TWV , T T ., 

FQRMATC20X, 'MIN C(I,J)*X(I,J) ' ,F10,4, 'CORRESPNG D(I, J)*X(I,J) ' 

1 ,5X,FIQ»4) 

RETURN 

WRITEdOUT, 10005) - #yyCA #lllw . , . 

FORMATC50X, 'The primal is in£easbible'/750x, 'The problem is abor 
2ted and next problem pursued' ) 

RETURN 

WRITE(IQUT,1105) 

RETURN 

FORMaI(§OxI U WHOLE RANGE COVERED') 

RETURN 

END 


SUBROUTINE FEACHKC INFCON) 

CQMMON^/Bl ,MA (50) ,MB(50) , ICOSTt (50,50) ,ICOST2(50,50) 

MSHRT=0?MSUP=0 

INFCON=,TRUE, 

DO 10 I * 1 ,M 
MSUP=MSUP + MAC I) 

MSHRTaO 

IF( (MA(I) IlT.O) .OR. (MB(J) *LT, 0) .OR. (ICOSTi (I,J)»LT,O),0R* 


MSHRT ........ - . 

IF(MSHRT.NE.MSUP) 
INFCON*, FALSE, 
RETURN 
END 


GO TO 15 


505 FORMA 


g!BB8B*JK/SS? ?»i apSWSS?! «-«o . 2 , . 
ls¥6g?/|W T ' l6uT 

LOGICAL ILEMTY 
ILEMTY», FALSE, , 

WRITE(lSUT,505) jILCNT 



15) 



48500 

48600 

48700 

48800 

48900 

49000 

49I0Q 

49200 

49300 

49400 

49500 

49600 

49700 

49800 

49900 

50000 

50100 

50200 

50300 

50400 

50500 

50600 

50700 

50800 

50900 

51000 

51100 

51200 

51300 


C 

c*** 

c 


DO 100 I s l, 

TF(IL(I,1).GT 

CONTINUE 

ILEMTX=,TRUE, 

RETURN 

IPQSTNsI 

RETURN 

END 


1,ILCNT 
GT»0) GO 


SUBROUTINE XCHANG CPSOL, COST »X , M » SUM, 

COMMON /B1 /NOSP , NOSHRT # MA (50 ) # MB ( 50 , 

DIMENSION PSOL(N) , COSTf N) ,X(5Q,5Q) 

INTEGER ROW, COL 

SUM=0,0 

DO 30 I = 1 , N 

IPxpSOL(I)+0,5 

PSOLCI )*IP 

RQW=(I-1)/N0SHRT +1 

COL-I" ( (ROW-1 ) *N05HRT ) 

SUM=SUM 4 PSOLC I ) *COST( I ) 
X(ROW,CQL)=PSOL(I) 

CONTINUE 

RETURN 

END 


,ICOST1(50,50), IC0ST2 (50 


B- 8 



APPENDIX C 


ES2SI32 isi S'BS£2£i2B 21 ESD222 gJSSiSBS 


00100 
00200 
00300 
00400 
00500 
go 600 
00700 
00800 
00900 
oiooo 
oiioo 
01200 
01300 
01400 
01500 
01600 
01700 
01800 
01900 
02000 
02100 
02200 
02300 
02400 
02500 
02600 
02700 
02800 
02900 
03000 
03100 
03200 
03300 
03400 
03500 
03600 
03700 
03800 
03900 
04000 
04100 



05000 

05100 

05200 

05300 

05400 

05500 

05600 

05700 

05800 

05900 

06000 

06100 

06200 

06300 

06400 

06500 


c 

c 

C 


4 were generated, 
by 5 were generated, 
by 3 were generated, 


C 

c 

c 

c 


10 


12 


15 


16 

200 


205 


1949 

18 


I 


31 


35 


40 


With ISEED 2948513, ten problems of 4 by 
wltn ISEED 4548511 , twenty problems of 5 
with ISEEp 294851J, twenty problems of 3 
IMPLICIT INTEGER (A-Z) 

REAL N i(5000) PPLY(50) ' DEMAND(50),COST,C50#50) ' COST2C50#50) 
NQQFPR=10 ; 

I0UT=21 

tcc ,_„OPEN CUN1T=I0UT,FILE='GENERT,0UT') 

f|||S given for the first set of problems of 4 by 4 is 2948513 
ifEED given for the second set of problems of 5 by 5 is 4548511 
ISEED GIVEN t OR THE FOURTH SET OF PROBLEMS OF 7 AY 7 IS 40205096 
TSEED FIFTH SET OFPROBLEMSOFe BY*6 I S 9 5 541 7 1® 

NOTIME=l 

DO 95 I = 1, NOTIME 
IFQ.EQ.USIZE = 6 

IFa.NE ,i5size=size* 2 
N=2*<SIZE**2)+2*SIZE 
DU 90 JJJslNOQFPR 
DO 10 II s I ,n 
R(II)=0,0 
ISEED3ISEED+1 
CALL GGUB(ISEED,N,R) 

TOTALsO? RINDEXsO 
DO 15 II = 1 .SIZE 
RINDEX=RINDE&t 1 
INTERsRCP INDEX)*! 00.0 

IS4LY?f!f£|SfEg* INTER - GT - ,5)G ° T ° 12 

C0NTl5S| TAL + SUPPLYai) 

PAUSE 'FINISHED SUPPLY' 

DO 20 II=»i , SIZE 
IFQIssSIZ^JGO TO 19 
RINDEXsRXNDEX+1 

lU¥o%&A 20a ' 20S 

&^ B(ISEED ' N ' B) 

GO TO 16 



Tl}T ALssTOT 

Jo < TO T 20 • E9<1,c * u ' m 0DIFY (DEMAND , II (TOTAL) 

S5MS8AJ T, " ,nm 

PAUSE 'Finished demand' 

DO 50 11=1 , SIZE 
DO 50 JJ=1 * SIZE 
DO 45 J = 1 ,2,1 
RINDEX=RINDEXtl 

||fe|El|S BCISEED ' N ' R> 

S»3:I8:U S8Bmi:SH:iSI6 



06600 

45 

06700 

06800 

50 

06900 

110 

07000 

07100 

07200 

07300 

07400 

07500 

07600 

101 

07700 

82 

07800 

c 

07900 

90 

08000 

08100 

08200 

08300 

08400 

08500 

08600 

08700 

08800 

08900 

95 

09000 

09100 

5 

09200 

09300 

09400 

09500 

09600 

09700 

10 


CONTINUE 

CONTINUE 

WRITE( I OUT #110) JJJ, SIZE, SIZE 

FORMAT! ' Problem no. ',14, ' of size ',13, ' x ',13) 
FQRMATC40I3) 

WRITE! I OUT, 101 ) ! SUPPLY (II), 11 = 1, SIZE) 

WRITE! IOUT, 101) (DEMAND (II), 11=1, SIZE) 

DO 82 JN=1 ,2 
DO 82 JP=1 , SIZE 

IF(JN==l)WRlTE!I0UT,101) ICOSTl !JP,II),II=1, STZE) 

IF ( JN==2 ) WRITE! IOUT, 101 )!COST2(JP, II), 11=1, SIZE) 
CONTINUE 

PAUSE ' One problem is over' 

CONTINUE 

CONTINUE 

STOP 

END 

SUBROUTINE MODIFYIDEMAND, II , TOTAL) 

IMPLICIT INTEGER !A-Z) 

DIMENSION DEMAND!50) 

MAX=DEMAND!1) 

DO 10 1=1,11 

IF ( DEMAND !I),GE.MAX) GO TO 5 

GO TO 10 

MAX=DEMAND!I) 

INDEX = I 

CONTINUE 

DEMAND ! INDEX) =1 

TEMP = MAX-DEMAND ! INDEX) 

TOTAL=TOTAL + TEMP 

RETURN 

END 


C- 2 
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